Multi-particle entanglement manipulation under positive partial transpose preserving 

operations 



Satoshi Ishizaka^'^ and Martin B. Plcnio^"'* 

^Fundamental and Environmental Research Laboratories, 
NEC Corporation, 34 Miyukigaoka, Tsukuba, 305-8501, Japan 
^PRESTO, Japan Science and Technology Agency, 4-1-8 Honcho Kawaguchi, 332-0012, Japan 
'^QOLS, Blackett Laboratory, Imperial College London, Prince Consort Road, London SW7 2BW, UK 
^Institute for Mathematical Sciences, Imperial College London, 53 Exhibition Road, London SW7 2BW, UK 

(Dated: February 1, 2008) 

We consider the transformation of multi-partite states in the single copy setting under positive- 
partial-transpose-preserving operations (PPT-operations) and obtain both qualitative and quantita- 
tive results. Firstly, for some pure state transformations that are impossible under local operations 
and classical communication (LOCC), we demonstrate that they become possible with a surpris- 
ingly large success probability under PPT-operations. Furthermore, we clarify the convertibility 
of arbitrary multipartite pure states under PPT-operations, and show that a drastic simplification 
in the classification of pure state entanglement occurs when the set of operations is switched from 
LOCC to PPT-operations. Indeed, the infinitely many types of LOCC-incomparable entanglement 
are reduced to only one type under the action of PPT-operations. This is a clear manifestation 
of the increased power afforded by the use of PPT-bound entanglement. In addition, we further 
enlarge the set of operations to clarify the effect of another type of bound entanglement, multipar- 
tite unlockable bound entanglement, and show that a further simplification occurs. As compared 
to pure states a more complicated situation emerges in the mixed state settings. While single copy 
distillation becomes possible under PPT-operations for some mixed states it remains impossible for 
other mixed states. 



I. INTRODUCTION 

Constraints and resources are intimately related in 
physics. If we impose a constraint on a physical set- 
ting then certain tasks become impossible. A resource 
must be made available to overcome the restrictions im- 
posed by the constraints. By definition such a resource 
cannot be created employing only the constrained set of 
operations but it may be manipulated and transformed 
under these operations. That the amount of resource 
does not increase under any operation satisfying the con- 
straint emerges then as a fundamental law, for example 
in entanglement theory 

One example of particular importance is the restriction 
to local quantum operations and classical communication 
(LOCC). The resource that is implied by this constraint 
are non-separable states and in particular pure entangled 
states such as singlet states, neither of which can be cre- 
ated by LOCC alone. This setting gives rise to a theory 
of entanglement as a resource under LOCC operations. 

Any such theory of entanglement as a resource will 
generally aim to provide mathematical structures to al- 
low answers to three questions, namely (1) the character- 
ization of entanglement, (2) the manipulation of entan- 
glement and (3) the quantification of the entanglement 
resource under the given constraint. Of particular 

interest is the question of how many incquivalent types 
of entanglement exist within such a theory. In the limit 
of infinitely many identically prepared copies of bipar- 
tite pure states, entanglement can be inter-converted re- 
versibly Q and it is reasonable to say that there is only 
one type of pure bipartite entanglement. Even for pure 



states, the situation changes dramatically when we con- 
sider the single copy setting. It has been shown that the 
Schmidt rank of bipartite pure states cannot be increased 
by LOCC Him 0,0. At the single copy level, the con- 
vertibility of bipartite entanglement is then characterized 
by the Schmidt-rank Q. For finite dimensional systems 
a state can be converted to another one with finite prob- 
ability exactly if the Schmidt-number of the target state 
is not larger than that of the initial state. In a tripartite 
setting the situation is more complicated. Here it is well- 
known for example that a GHZ state cannot be trans- 
formed to a.W state and vice versa 0. These states are 
then said to be incomparable. It can be shown that there 
are two incomparable types of tripartite entanglement in 
three qubits systems. The situation is even more compli- 
cated in multipartite settings composed by many parties 
[lOl IllL lia. [la . or infinite dimensional bipartite sys- 
tems |l5lll6j. where there are many (possibly infinitely 
many) incomparable types of entanglement. 

A different setting is presented by the concept of par- 
tial time reversal or partial transposition |l7j| . For two 
qubits, states that remain positive under partial transpo- 
sition (denoted as PPT-states) are exactly the separable 
states [l^ but for higher dimensions this is generally not 
the case as there are PPT-states that are inseparable . 
This motivates the definition of the set of positive-partial- 
transpose-preserving operations (PPT-operations), de- 
fined as operations that map any PPT-state into another 
PPT-state 20i. In this case, the resource are states that 
are not PPT (denoted as NPT-states). In the single copy 
setting for pure states, it has been shown that both un- 
der PPT-operations l2ll| and with LOCC supported by 
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PPT-bound entanglement the Schmidt-number can 
be increased so that state transformations become pos- 
sible that are strictly impossible under LOCC. Further- 
more, there are mixed state transformations that are re- 
versible in the asymptotic setting . This suggests that 
a theory of entanglement under PPT-operations might 
have a much simpler structure than that under the LOCC 
constraint. 

In this paper, we focus attention on the entangle- 
ment manipulation under PPT operations in the non- 
asymptotic, single copy setting to explore what simpli- 
fication occur. We consider PPT state transformation 
in multipartite settings and obtain both qualitative and 
quantitative results. In Sec. the general settings 
and notations of PPT preserving operations are intro- 
duced. In Sees. mil and llVI we first demonstrate that the 
transformations of pure states that are impossible un- 
der LOCC become possible with a surprisingly large suc- 
cess probability when employing trace preserving PPT- 
operations. In Secs.E a rather tractable scheme of trace 
non-preserving PPT-operations is introduced and dis- 
cussed. We will then completely clarify the convertibility 
of all multipartite pure states under PPT-operations in 
Sec. IVII In Sec. IVIII we enlarge the set of operations 
beyond that of PPT-operations to consider the effect of 
multipartite unlockable bound entangled states. In Sec. 
IVIIII we will consider the transformation of a single copy 
of mixed states into pure entangled states, i.e. the single 
copy distillation under PPT-operations. A summary and 
conclusion is given in Sec. IIXI 

II. BASIC NOTATION 

To begin with, let us denote Ti-iV) {H{V')) the space 
of Hermitian operators on the Hilbert space V {V). A 
superoperator ^> from V to V' is a linear transformation 
from TCiV) to There is a natural isomorphism .2(1] 

which associates with superoperators 4' : ?i{V') 
a Hermitian operator ri(^') e HiV) (g) H{V') such that 
for all A e niV) and B e niV) we have 

tr{*(A)B} = tr{rj(*)yl » B}. (1) 

Maps that are trace non-increasing then satisfy 

tiV'M'^)} < ly (2) 

with equality if 5* is trace preserving. A superoperator 
'i' is called positive if for any A > we have ^'(A) > 
and it is called completely positive if IvF ® ^ > for any 
space W. Following j23| complete positivity (CP) of ^' 
can be verified by checking 

> (3) 

where Ty denotes the partial transposition with respect 
to V. 

An additional concept comes into play when we con- 
sider multipartite systems. A CP-map on bipartite sys- 
tems * : niVA) ® niVB) -> niVX) ® niVl^) is called 



positive partial transpose preserving (PPT) |23|, if we 
have Pa o vl/ o > (Pb o * o P^ > 0) for the partial 
transposition map Ta (Pb) with respect to party A (B). 
On the level of the state f^('I'), this condition reads 

(17(^')r^)^^-*®^^A > or (f7(1')^^)^^«®^^A > 

where Py^ (Py ) denotes partial transposition applied 
to space Va (V^)- In the bipartite case, there are two 
equivalent choices for the partial transposition. In the 
tripartite setting however, there are three different possi- 
ble partial transpositions that are generally not equiv- 
alent. A CP map ^ : 7^(Va) ® H(Fb) «) W(Vb) ^ 
n{VX) ® H(F^) ® n{Vl;) will be called PPT in the fol- 
lowing if 

(f7(*)^^)^^-®^n' > (4) 

for all i — A,B and C. 

Let us now consider the transformation of a state 
p e T^lV) into a state a E H{V') with the probabil- 
ity of p{p — > a). For this probabilistic transformation, 
we construct the trace preserving CP-PPT map with 
two outcomes, one that gives a and one that gives some 
other state. The two parts are given by CP-PPT maps 
^' and ijj, respectively. The associated Hermitian opera- 
tors are denoted by fl and uj. The map ^' then satisfies 
"^{p)—p{p^a)a or 

tr{1'(p)(l - a)} = ti{n{^)p ® (1 - cr)} = (5) 

when tr is a pure state. The success probability is then 
given by 

p(p ^ cr) = tr{«'(/9)} = tr{ri(^')p (g) 1} 

tr{f2(*)/9(K)o-}. 

The PPT-map ijj, on the other hand, does not suffer any 
constraint other than the condition of trace-preservation 
for "if + ip. On the level of states, the trace-preserving 
condition is 

trv + uj} = Iv , (6) 

where, as we will do in the remainder of this paper, we 
have dropped the ^E" {tp) in il(^') {uj(ip)) for brevity. It 
should be noted that a rather simple structure can be 
assumed for w without loss of generality. Let us con- 
sider a map x which maps arbitrary states in HiV) into 
a maximally mixed state of ly'/dim{7i(l/')} £ ^{{V'). 
This map is a trace-preserving CP-PPT map since the 
corresponding state is ly ® l\/'/dim{H(y) H{V')}. 
Therefore, a composed map of x o i/; is a CP-PPT map 
if t/; is a CP-PPT map. Furthermore, if ^ + is trace 
preserving, ^f + x o V' is also trace preserving, and hence 
the replacement of "0 by x o ^/j does not alter p{p^ <j). 
One may then assume ■(/) = x o V' since the output of tjj is 
arbitrary. On the level of the state, this assumption is 

" = "^^ dim{HV)} - 
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In the subsequent Secs. lIIII and llVI we maximize p{p^a) 
for some important classes of pure states in both bipartite 
and tripartite settings. In particular, we demonstrate 
that transformations of pure states that are impossible 
under LOCC can be achieved under PPT operations with 
a surprisingly large success probability. 



III. CONVERSION OF MAXIMALLY 
ENTANGLED STATES 

For two d-dimensional systems we denote the maxi- 
mally entangled state by = I^J I where 

d-l 

^ 73 SI'"* 



i=0 



In the single copy setting, it is known that LOCC cannot 
increase the Schmidt rank of a pure state 0, 0, S S 13 ■ 
Therefore, p{Pj^ P^) = for LOCC transformation 
whenever d! >d. 

In the following we proceed with the construction of 
the CP-PPT maps 5" and -0 that maximize the suc- 
cess probability for this transformation. For d! > d this 
amounts to the maximization of 



P{Pl ^ Pt') = ^^{^Pl ® Pi'} 
under the constraints 

ii{Q.P+ ® {t - P+)} = 0, try^O + cj} = 1, 



(8) 



> 0, 



if^ > 0, (9) 

CJ^^ > 0, 



where P+ eU{V) and P+en{V'). Since both Pl^Pj, 
and — Pj) are invariant under the local unitary 

transformation of Ui ® Ul ® U2 ® U2 with Ui and U2 
being arbitrary unitary operators, it suffices to consider 
VL and Lo that are invariant under these local operations, 
i.e. 

^ = aiP+ ®P+ + 02(1 -P+)®P+ 

1 — P+ 1 — P+ 

+ "sP/ ® + 04(1 - P/) — 
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The first two constraints in Eq. © yield 03 = 0, &i = 
1 — oi, and 62 = 1 — a2 — 04. These equalities can be used 
to eliminate bi and &2 in the remaining constraints. The 
remaining constraints then result in 

1 > Ol > 0, 02 > 0, 04 > 0, 1 > 02 -I- 04, 

(d' + l)ai + (rf' + 1)02 + 1)04 > 0, 

-(d' + l)ai + (rf' + 1)02 + 1)04 > 0, 

-(d' - l)ai - (d' - 1)02 + (d- 1)04 > 0, 

(d' - l)ai - (d' - l)(d+ 1)02 + (d+ 1)04 > 0, 

-ai - (d - 1)02 - (d - 1)04 + d > 0, 

ai - (d+ 1)02 - (d+ 1)04 + d > 0. 



The constraints in the first row arise from > 
and > 0. The last two rows are due to 

(^^'">'a®'"v^ > Q g^jj^j ^lyQ remaining for inequalities 
arise from (^^^^j^'^^^a'^^v > ^j^g maximization of 
p(P^—>P^)=ai under these constraints is a linear pro- 
gram and we can identify the optimal solution as ai = 
d(d-l)/(dd'-M'-2d), 02 = 0, and 04 = d(d'-l)/(dd'-M'-2d). 
Consequently, for d' > d the optimal probability for the 
transformation of P^ into P^, thereby increasing the 
Schmidt rank, under PPT-operations is given by 



PiPl - Pi) 



d(d- 1) 
dd' + d' - 2d' 



(10) 



We emphasize that this success probability is nonzero 
even when d' > d > 2, while it is strictly zero for the 
LOCC transformation. 



IV. CONVERSION FROM GHZ TO W STATE 

In the tripartite setting, it is well-known that the suc- 
cess probability p{GHZ—>W) = for the LOCC trans- 
formation from a single copy of 



\GHZ) = 



1 000) 



1111) 



to 



V2 

|001) -f- 1010) + |100) 
V3 



[9|. In the following we will demonstrate that this is not 
the case when we consider PPT-operations. To this end, 
we maximize 

P{PGHZ Pw) = tl{Q,pGHZ ® Pw} (11) 

under the constraints for i — A, B, C, 

tr{npGHZ (^{t- pw)} = 0, trv'{n + Lo} = 1, 
(•^rv,»r;,^)rv > g, ff^^ > 0, 



w > 0, 



where poHZ = \GHZ) {GHZ\ eH{V) and pw = \W){W\E 

n{v'). 

The solution of the problem is greatly aided by the 
use of a number of symmetries. Indeed, both the states 
Pghz®{^ — Pw) and pghz®Pw are invariant under the 
local operations 

(a) X®X®X®t®t®t, 

(b) Z(g)Z(g)l(g)l(E)l(g)l, 

(c) IL(8)Z(g)Z(8)]l(g)]l(8)l, 

(d) l(g)l(E)l(g)Z(g)Z(g)Z, 

(e) Pi Pi Pi 1 (g) 1 a , 

(/) a «) 1 (g) a «) P2 «) P2 «) P2, 
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where Pi = |0)(0| + |l)(l|e2"/3 and P2 = e"/2|0)(0| + 
|l)(l|e'^\ These local symmetries are supplemented by 
the non-local joint permutation symmetry 

(g) ^(123) X ^(456) 

where V represents an arbitrary index permutation. The 
symmetries (a) - (g) allow for a considerable simplifica- 
tion of n and oj. Indeed, the symmetries (b), (c) and (e) 
ensure that the matrix elements ^ij^jj^kihmim,i232k2i2m2n2 
can only be non-zero if the indices satisfy simultaneously 
h =^2, ji =32 and ki =k2 or ii = l-?2, ji = l-j2 and 
ki = l — k2. The symmetry (g) yields 



abcdef.ghijkl 



V{abc)V(deS),V{ghi)V{jkl) 



(12) 



for any index permutation V . Symmetry (a) yields 

^000/imini ,000/2"^2n2 ~ ^lll/lmlnl,lll/2^^^2^^2 5 (1^) 
^001/imini ,001/2Tn2n2 ~ ^110/i7ni ni ,110/2m2n2 ) (-^^) 
^000/imini,lli;2m2n2 ~ ^lllhmi "l .000/2m2n2 • (l^) 

Presenting all nonzero matrix elements of ^abcdef,ghijki 
for {abc,ghi) = (000,000), {abc,ghi) = (001,001) and 
{abc, ghi) = (000, 111) fixes all other matrix elements by 
virtue of the symmetries Eqs. (|12I15I) and the Hcrmiticity 
oi n. To obtain a trial solution we chose 

^000000,000000 = ^001000,001000 = ^^000000,111000 = bi, 

^000001,000001 = f^OOOOOl, 000010 = ^000001,000100 = ^2, 

^^000011, 000011 = ^^000011, 000101 = ^^000011, 000110 = ^4, 

f^OOlOOl, 001001 = —f^OOlOOl, 001010 = —f^OOlOOl, 001100 = &2, 

f^OOlOlO, 001010 = f^OOlOlO, 001100 = f^OOllOO, 001100 = b2, 

f^OOlOll, 001011 = f^OOlOll, 001101 = —f^OOlOll, 001110 = — &4, 

f^OOllOl, 001101 = —f^OOllOl, 001110 = f^OOlllO, 001110 = bi, 

f^ooiiii, 001111 = Siloooiii, 000111 = — 3rioooiii,iiiiii = Sfeg, 

f^OOOOOl, 111001 = f^OOOOOl, 111010 = f^OOOOOl, 111100 = ^2, 

f^OOOOlO, 111010 = f^OOOOlO, 111100 = f^OOOlOO, 111100 = ^2, 

f^OOOOll, 111011 = f^OOOOll, 111101 = f^OOOOll, 111110 = —bi, 

f^OOOlOl, 111101 = f^OOOlOl, 111110 = f^OOOllO, 111110 = —bi. 

Likewise, the non-zero matrix elements of Luy can be con- 
structed from 

(t^y)ooo,ooo = 1 — fee — 3^4 — 3^2 — ^1, 
(t^y)ooi,ooi = ('^v)ooo,iii, 
(t^y)ooo,iii = 66 + 364-362 + 61, 

where we chose 

61 = 



1 + Vl - 4x2 ^ X ^ 62 9^2 

^ , 62 = ^, 64 — — , be — — , 

D 3 Oi 3x 



X = i(-2 + (18-6\/3)i/3 + (18 + 6\/3)i/3). 
8 

A lengthy but elementary calculation (preferably exe- 
cuted employing a program capable of symbolic manipu- 
lations) then confirms that this trial solution satisfies all 



the constraints and yields the success probability 

ti-{npGHZ <E) pw} = 662. (16) 

We then consider the dual problem of the primal prob- 
lem Eq. ((TT|l psj l- Every feasible point of the dual prob- 
lem provides an upper bound on the solution of the pri- 
mal problem Eq. pi|) . The above result of Eq. (|16|l is 
then proven to be optimal as shown in Appendix 1X1 

As a consequence, the optimal probability for the 
transformation of a GHZ to a state under PPT- 
operations is given by 



p{GHZ -^W)^ 662 w 0.75436. 



(17) 



that is more than 75%. This very high success probability 
is somewhat surprising, since the success probability for 
the LOCC transformation is strictly zero. Note that this 
result also implies that a GHZ state can be transformed 
into a W state employing LOCC supplemented by PPT- 
bound entanglement. 



V. TRACE NON-PRESERVING CP-PPT MAPS 

In the previous two sections we have demonstrated 
explicitly that the success probability for the transfor- 
mation between pure states can in some cases be im- 
proved significantly by employing PPT-operations in- 
stead of LOCC operations. Obtaining the optimal suc- 
cess probabilities is a hard task, however, especially in 
the multipartite setting. In the following we will consider 
the slightly more tractable setting of trace non-preserving 
PPT maps. In this setting we also optimize a CP-PPT 
map 5* or equivalently the associated state fi, but the 
trace preserving condition of Eq. ^ is replaced by the 
trace non-increasing condition of 



V- 



(18) 



As a result, the completion ip of the map is a CP map 
but it is not necessarily a PPT map. This will gener- 
ally allow to find success probabilities for state trans- 
formations that are larger than those obtained under 
trace-preserving PPT operations. It is important to note, 
however, that any transformation that possesses a non- 
vanishing success probability under trace non-preserving 
CP-PPT maps will also have a non-vanishing success 
probability under trace preserving CP-PPT maps. To 
see this, let il(^) be the state corresponding to a trace 
non-preserving CP-PPT map 5*. Since the completion 
-0 is not necessarily a PPT map, uj{^)^^ is sometimes a 
NPT-state. However, if we consider the states of i7'(^'') = 
en{^) and uj'{ij')^euj{'ilj) + {l-e)t(E)l/dim{niV')}, the 
state {oj')^^ becomes a PPT-state for a nonzero value of 
1 > e > 0. Both (f]')^'' and (w')^'' are PPT states sat- 
isfying the trace preserving condition of Eq. lO, and 
accomplishes the same transformation as albeit with a 
smaller success probability. In this way, one can always 
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construct a trace preserving CP-PPT map from the trace 
non-preserving CP-PPT map giving the same transfor- 
mation. 

The optimal probability in the trace non-preserving 
scheme for the transformation of maximally entangled 
states {d! > d) can be obtained in the same fashion as 
section III. Employing the notation of section III we ob- 
tain the constraints 



1 > ai > 
{d' + l)ai - 
-{d' + l)ai + {d' + l){d 



(d'-hl)(d-l)a2 + (d-l)a4 > 0, 

-id' + l)ai + id' + l){d + l)a2 + (d + l)ai > 0, 

-{d' - l)ai - [d' - l){d - l)a2 + {d- 1)04 > 0, 

(d' - l)ai - (d' - 1)02 + (d+ 1)04 > 

under which the success probability, given by ai, has to 
be maximized. The result is 



d-i 



d'~V 

whose PPT map ^' is, on the level of the state fl, 
d-1 1 



(19) 



n = 



Y-r\Pd®PI'+^j2^i^-PI)®i^-PI')- (20) 



It is noteworthy that the probability of Eq. (|19|l can be 
written as a ratio of the negativity of the initial and target 
state, i.e. 

where N{a) = (tr|cr^| - l)/2 [H This somewhat 

fascinating expression resembles the case of the LOCC 
transformation of pure states, where the optimal proba- 
bility agrees with a ratio of a LOCC-monotone such that 
the partial summation of squared Schmidt coefficients @ . 
Although the monotonicity of the negativity in trace non- 
preserving PPT-operations has not been proved yet (in 
trace preserving PPT operations with a single outcome 
the negativity is a monotone llj), the tractable expres- 
sion of Eq. (|19|l is likely to be explained as a ratio of some 
monotone function. 

In the tripartite setting, the optimization of the success 
probability is still a hard task even in this trace non- 
preserving scheme. The result of the optimization for 
the transformation of GHZ — > is 



p{GHZ ^ W) 



5' 



and for the transformation of — > GHZ we have 
p{W ^ GHZ) = i 



(21) 



(22) 



The proof for these two results are described in appen- 
dices El a-nd [ni This result implies that the transforma- 
tion of ^ GHZ is also possible by trace preserving 
PPT-operations, although the optimal probability may 



be smaller than 1/3. Therefore, PPT-operations can 
inter-convert even the LOCC-incomparable pure states. 
In the next section, we will completely clarify the con- 
vertibility by PPT-operations for all multipartite pure 
states in the single copy setting. 



VI. CONVERTIBILITY OF PURE STATES 

In this section we will consider the transformation be- 
tween single copies of A^-partite pure states under PPT- 
operations. By definition, PPT-operations map PPT- 
states to PPT-states. As a consequence, transforma- 
tions such as \(I)ab)'^\^c) ~^ \GHZ) or \(j)\g) '^\Qc) ^ 
\Qa) ® l^ijc) ^'"^ impossible, since they are not PPT- 
preserving with respect to party C. Therefore, let us 
first assume for the transformation of | '0) — » 1 0) that both 
IV') and are "genuinely" entangled over all A^ parties. 
This assumption means that 



mm^'t^ and (|0)((/.|)r» ^0, 



(23) 



for all possible bipartite partitioning of i. For exam- 
ple, i = A, i?, C in a tripartite setting, and i = A, B, 
C, D, AB, AC, AD in a four-partite setting. As dis- 
cussed in the previous section, it suffices to consider trace 
non-preserving CP-PPT maps ^' in order to check the 
convertibility under trace preserving PPT-operations. 
Therefore, we will construct an f2 satisfying the con- 
straints 



tr{r!|^)(V'|®(l-|0)(0|)}-O, 



(24) 



where £T-C{V), |0) ETi.{V'), and i stands for any pos- 
sible bipartite partitioning as explained below Eq. H23|l . 
We have omitted the trace non-increasing condition, be- 
cause we are not interested in the explicit value of the 
success probability but only whether it is zero or not. In 
view of Eq. H2UI) , a suitable trial form is 

n = ® |0)(0I + (1 - WW) ® (1 - I'/'X'^I), (25) 

for which the first two constraints in Eq. 1)24(1 are satisfied 
when x>0. Furthermore, due to the assumption of Eq. 
(123), the last constraint {fl^''^'^'')^^ > is also satisfied 
for an appropriate value of a; = xq > as shown in |22j | . 
As a result, for x = xo we have 



tr{r!(|^)(V'|® 10) (</>!)} =a:o >0, 



(26) 



so that for arbitrary pairs of genuine A^-partite entan- 
gled states of IV') and we can always find an il such 
that p(IV') ^ 10)) > 0. As a consequence, all genuine 
A^-partite pure entangled states are inter-convertible by 
PPT-operations. In this way, the classification of A^- 
partite entanglement is drastically simplified when we 
consider PPT-operations. 

Let us next investigate the convertibility between a A^- 
partite state IV'^^^) and a (A^- l)-partite state |0(^-^)). 
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iV-partite 



Tripartite 

(A'=3) 

Bipartite 

{N=2) 



Separable 





PPT-operations ■ 
unlocl<able BE < 



can conclude that the transformation such as GHZ ^ W 
can be accomphshed by LOCC with the consumption of 
PPT bound entangled states. Much attention has been 
paid to bound entanglement to clarify its properties, and 
several applications of bound entanglement have been re- 
ported m El US El El El Ei El El El El El 

As shown above, PPT bound entanglement enables the 
LOCC implementation of large classes of entanglement 
transformations that are impossible by LOCC alone. 



VII. UNLOCKABLE STATES AND 
CONVERSION OF PURE STATES 



FIG. 1: The classification and convertibility of multipartite 
pure entangled states under PPT-operations. r denotes the 
Schmidt rank of bipartite entanglement, and the set of entan- 
gled parties in (A'^— l)-partite entanglement is assumed to be 
a subset of the set of entangled parties in A^-partite entan- 
glement. There is only one type of A'"-partite entanglement 
under PPT-operations. Furthermore, the convertibility with 
the support of unlockable bound entanglement (BE) is also 
shown (see also ^3)- 

It is obvious that l^'-^^"'^-') — > IV'^^'') is impossible be- 
cause such a transformation is not PPT-preserving. Like- 
wise the transformation of |?/)(^)) — > \(f>^^~^'>) is impos- 
sible if the set of entangled parties in \(f>^^^^'>) is not 
a subset of the set of entangled parties in (e.g. 

\''Pabc) ^\'^ad) impossible). Otherwise, the transfor- 
mation is possible because a A^-partite GHZ state can be 
transformed to a (A^— l)-partite GHZ state by LOCC, 
and hence the sequential transformation of IV'^^'') ~* 
IGiTZW) ^IGilZ(^-i)) ^ is possible. The 

classification and convertibility of arbitrary multipartite 
pure entangled states under PPT-operations are summa- 
rized in Fig. ^ 

It is important to note here that the power of PPT- 
operations, by which A^-partite pure entangled states be- 
comes inter-convertible as discussed above immediately 
implies that the same holds for LOCC supported by PPT 
bound entanglement. This is due to the fact that any 
PPT-transformation can be accomplished (with smaller 
but nonzero probability) by LOCC supported by the ad- 
ditional resource of PPT-states j25j (see the note of [2^). 
Indeed, 

*(p) = trvM'i'f'' ip^'' ® 1)}, (27) 

and the state of fi(^')^^ >0, which is a PPT-state if "if is 
a CP-PPT map due to {n{^f^^^'')^''>0, is utilized and 
consumed in the LOCC-implementation of ^'(p) ^2511. If 
a CP-PPT map "if can accomplish a transformation that 
is impossible under LOCC alone, then ri(^)'"^ must be 
entangled (otherwise the transformation can also be ac- 
complished by LOCC because LOCC can generate any 
separable state), and therefore the PPT-state 17(5')'"^ is 
a PPT bound entangled state Consequently, one 



As mentioned in the previous section, the transforma- 
tion 

I'/'ab) ® I0c) ^ \GHZabc) (28) 

cannot be achieved even when PPT-operations are em- 
ployed, and therefore cannot be achieved by LOCC sup- 
ported by PPT-bound entanglement. However, it has 
been shown that a GHZ state can be distilled from a 
tripartite NPT-bound entangled state, if A and B per- 
form a global operation on the state [Hi. Such NPT- 
bound entangled states are called unlockable states be- 
cause bound en tang lement is unlocked by the global op- 
eration m m Ef. The global operation of A and B 
can be accomplished by LOCC consuming \4>\g), and 
consequently the transformation of Eq. (|28|l is possible 
when LOCC is supported by the unlockable bound en- 
tanglement [111. Likewise, unlockable states which can 
be utilized for the LOCC-transformation from a (A^— 1)- 
partite GHZ state to a A^-partite GHZ state have been 
shown in [29l |. In this section we consider this type of 
transformation using a certain general scheme. 

To this end, we generalize PPT-operations by relax- 
ing the PPT-preserving condition with respect to C, 
^^rvc^r^^^jFv >o^ which is responsible for the impossi- 
bility of the transformation of Eq. (|28() . We will therefore 
construct a map ^ whose associated Hermitian operator 
il satisfies 

ty{nP+g ® (1 - PGHz)} - 0, r!^^' > 0, 
(f^rv^^rv^^rv > (n^^s'^'^v^fv > (29) 

whereP+s-|0+s)(0+s|eH(T/) and pghz (^n{V'). As 
a trial form for fi, we adopt again Eq. H25|) . i.e. 

n = xPj,B ® PGHZ + (1 - Pab) ® (1 - PGHz)- (30) 

As mentioned in the previous section, 

{Pisf^^O and p^ahz^O 

ensure the existence of > such that (fl^^-^^^^A > 
for < X < xq (indeed, we have = 3), and like- 
wise with respect to B. We can now check easily that 
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all constraints in Eq. H29|l are satisfied for a; = 3, yield- 
ing a nonzero success probability p{\(t) 



ABI 



■\GHZ))>0 

since iv{Q.{P'^g ® Pghz)} > 0. Consequently, the trans- 
formation of \4'\g) \GHZ) is possible under the set 
of operations that maps NPT-BE states with respect to 
party C into itself, as expected. Employing symmetries 
'^^ ^AB ® PGHZ, the optimized success probability in the 
trace non-preserving scheme is then obtained as 



p(|0+5)^|Gi7Z)) = -, 



(31) 



and, on the level of states, the map 4* realizing this suc- 
cess probability is given by 



-p 



AB 



' Pghz 



:(1 



Pab) ® (1 



PGHZ~Pmi - Piio), (32) 



state P^i , i.e. the single copy distillation from a mixed 
state employing PPT-operations. 

Let us consider the antisymmetric Werner state which 
is defined as 



pa 



(34) 



J>1 



where P^ is the projector onto the antisymmetric sub- 
space of (g) C^, and = / For the 
transformation of cr^ — > P^, we can construct CP-PPT 
maps of 5* and its CP-PPT completion ij} employing the 
twirling symmetries of the two states. The result of the 
optimization is, on the level of the state (the state w 
is given by o^y = 1 — try f2). 



dd' +d' -2d 



Pa 



Pi' 



id' - l)P^ 



d'^ - 1 



where pooi = |001) (001| and pno = |110)(110|. It can be 
confirmed that the state O'"^ is an unlockable state as 
follows. Due to the constraints of {if''A<^'^W)^v^Q ^nd 
(yfvs^r'^^g-^Tv > 0^ the mixed state of n^"-' is undistill- 
able by LOCC, because LOCC is PPT-preserving and no 
tripartite and bipartite pure entangled state exists that 
is PPT with respect to both A and B. However, a GHZ 
state can be distilled from n^"" of Eq. ^ or Eq. if 
A and B perform global operations that distinguish P^g 

and l-P^B- 

Similarly, the map ^P, whose associated state is 



for d'>d>2, and 



d{d' ~ 1) 
2(d' + l) 
{d+l)d' 



pa I jd-d') p. 



ps 



{d + l)d' 



1 



Pi 



for 2<d' <d where PJ is the projector onto the symmet- 
ric subspace of (S) C^. The optimal success probability 
under trace preserving CP-PPT-operations is then given 

by 



^ = ^Pg'hz ' 



{N') 

' Pghz 



i^-p'Gflz)®i^-p'GWz), (33) 



(AT') 



for d'>d>2. 



can transform a A^-partite GHZ state (Pqhz) t*-* ^ 
partite GHZ state, and furthermore the state O'"^ is an 
unlockable state if TV' > iV > 2 ■ As shown in the pre- 
vious section, all genuine iV-partite entangled states are 
inter-converted by PPT maps. The composition of the 
PPT maps and the map given in Eq. (|33|1 is again a map 
whose associated state is an unlockable state. This im- 
plies that all pure entangled states can be inter-converted 
independently of the number of parties {N) when a single 
copy of an appropriate unlockable bound entangled state 
is available as a resource. In this way, the consumption 
of unlockable bound entanglement allows to overcome 
the LOCC-constraint between pure states with different 
sets of entangled parties, while the consumption of PPT 
bound entanglement overcomes the LOCC-constraint be- 
tween pure states with the same set of entangled parties 
(Fig.©. 



VIII. SINGLE COPY DISTILLATION 

So far, we have concentrated our attention on the dis- 
cussion of transformations between pure states. In this 
section, we will now consider the transformation of a sin- 
gle copy of a mixed state p into a maximally entangled 



pK^P1)-{ (35) 

Therefore, the success probability is nonzero for d'>2. 

On the other hand, the success probability for the same 
transformation under LOCC operations alone is strictly 
zero whenever d' > 2. This can be proven as follows: 
The I'tpij) in Eq. (|34|l are maximally entangled states on 
(8) C^. Therefore, each can be prepared from P^ 
by local unitary transformations only. As is an equal 
mixture of all possible \ip~j), cr^ can be prepared from a 
single copy of P^ by LOCC, and hence the transforma- 
tion of P^ has a finite success probability. If we 
furthermore assume that for d' > 2 the transformation 
a'^ — > P^ has a finite success probability under LOCC, 



P^ also has a finite 



then this implies that P2 cr^ 
success probability under LOCC. This contradicts that 
the Schmidt rank cannot be increased by LOCC. There- 
fore, the result of Eq. (|35|) implies that the success prob- 
ability of the single copy distillation is also significantly 
improved when PPT-operations are considered. 

It should be noted that the transformation of P2 
is possible under LOCC. Indeed, the local projection 

-|- |1)(1|, can accom- 
P^i is possible under 



P( 



P to crJJ, where P 



plish this. Furthermore, P2 
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PPT-operations, which enables the sequential transfor- 
mation of cr^ — * — > ■ Therefore, the feasibility 
of p{cr2 P^) can be regarded as being a consequence 
of the feasibility of p{P^ — > P^) under PPT-operations. 
Note however, that Eqs. 1(111}) and P5|) for d' > 2 imply 
that we have 

pK^P^)>p{a^,^P+)p{P+^P+). (36) 

Hence the direct transformation is accomplished with a 
higher success probability than that for the correspond- 
ing sequential transformation. 

The discussion above demonstrates that PPT- 
operations can improve the success probability of the 
single copy distillation for some mixed states. One may 
perhaps expect that single copy distillation becomes pos- 
sible for all NPT mixed states when we consider PPT- 
operations. This, however, is not the case. As shown in 
[4^ (see also 01 )i LOCC cannot distill any pure entan- 
gled state from a single copy of mixed states p on C^^C^ 
if rank(p) > cP—2. For such high rank mixed states, PPT- 
operations cannot distill any pure entangled state either. 
The proof of this statement is given in Appendix IdI 

This highlights the fact that LOCC state manipulation 
suffers certain restrictions that PPT-operations cannot 
relax. Indeed, the convertibility of some mixed states 
(into pure entangled states) at the single copy level, and 
therefore the convertibility of mixed states under PPT- 
operations remains much more involved than the convert- 
ibility of pure states. 



IX. SUMMARY 

In this paper we have considered the transformation 
of single copies of multi-particle entanglement under sets 
of operations that are larger than the class of local oper- 
ations and classical communication (LOCC). In partic- 
ular, we considered probabilistic state transformations 
under positive partial transpose preserving maps (PPT- 
maps). We demonstrated that transformations that are 
strictly impossible under LOCC can have a finite success 
probability under trace preserving PPT-maps. For spe- 
cific examples the optimal success probabilities are deter- 
mined. Surprisingly large values are obtained for example 
for the transformation from the GHZ to W state which 
under trace preserving PPT-maps has a success prob- 
ability of more than 75% while it is strictly forbidden 
under LOCC. Furthermore, we completely clarified the 
convertibility of arbitrary multipartite pure states under 
PPT-operations. As a remarkable result, we showed that 
all TV-partite pure entangled states are inter-convertible 
under PPT-operations at the single copy level, and there- 
fore infinitely many different types of entanglement under 
LOCC are merged into only one type. In this way, a dras- 
tic simplification in the classification of pure state entan- 
glement occurs when the constrained set of operations 
is changed from LOCC to PPT-operations. It should 



be emphasized that despite such drastic simplification 
in the single copy settings, the theory of entanglement 
under PPT-operations possesses the desirable properties 
that PPT-operations alone cannot create pure state en- 
tanglement and that the amount of bipartite pure state 
entan glem ent is uniquely determined in asymptotic set- 
tings 

The above results can be regarded as an application 
of PPT-bomid entanglement. In multipartite settings 
however another type of bound entanglement called un- 
lockable bound entanglement exists. Motivated by this, 
we enlarged the class of PPT-operations to consider the 
effects of unlockable bound entanglement. As a result 
we showed that all pure entangled states become inter- 
convertible independent of the number of parties, and 
therefore a further drastic simplification in the classifica- 
tion of pure states occurs when LOCC is supported by 
unlockable bound entanglement. 

Finally, we considered one aspect of mixed state entan- 
glement transformations, namely the single copy distil- 
lation by PPT-operations. We demonstrated that PPT- 
operations can distill a pure entangled state from a single 
copy of some mixed states with finite success probabil- 
ity, while the success probability under LOCC is strictly 
zero. However, we also proved that PPT-operations can- 
not distill pure entangled state from mixed states with 
very high rank. Therefore, certain restrictions of en- 
tanglement manipulation of mixed states under LOCC 
persist under PPT-maps, and the classification of mixed 
states under PPT-operations in the single copy settings 
is not as simple as that in the pure state case. 

It is important to further clarify how the structure 
of theory of entanglement is simplified under PPT- 
operations especially in the mixed state settings and in 
asymptotic settings, as this might enable a unified and 
systematic understanding of characteristics of quantum 
entanglement as a resource. 
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APPENDIX A: OPTIMALITY OF THE 
CONVERSION FROM GHZ TO W STATE 

In this appendix, we prove the optimality of Eq. H16|l . 
the probability for the transformation from GHZ to W 
state. To this end, we consider the dual problem of the 
primal problem Eq. I|ll|) '23*1 . The Lagrange function for 
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the minimization problem in Eq. is given by 

i=AM,C 

+tr{Ae(try/ {n + uj)-t)} + vtriflpGHZ «> (1 - pw)} 
+tT{\p{tTv'fl - 1)} + tr{Aep(trv"W - 1)}, 

where \,Xep,\A,)yBAc,tJ'A,lJ'B,lJ'C > 0. This La- 
grange function has to be minimized over all f2, a; > 0. 
This is feasible only if 

> ^ Af' + (pGHZ-Ap-Ae) (g) l-l^pGHZ {t-pw), 

i=A,B,C 

> -(Ae + Xep) ® 1 + A^A^ + pI"" + pIF, 

in which case we obtain the dual function 

g{Xp, Xep, Ae, ly) = -tr{Aep + Ae + Ap}. (Al) 

Every feasible point of the dual problem provides an up- 
per bound on the solution of the primal problem Eq. 
With the symmetries shown in Sec. II VI the Lagrange dual 
problem of primal problem Eq. is 

min tr{Aep -|- Ae -f Ap} (A2) 

under the constraints 

(PGHZ-Ap-Ae)(8)l-i^PGffz«'(l-PH')+ Y Af'''<0, 

i=A,B,C 

Xa,Xb, XcPa, /^s, A*c, Ap, Aep > 0, 
-(Ae + Aep) ® 1 P^/^ + Ms '"^ + Mc < 0- 

To prove the optimality of Eq. (|16|l . it suffices to provide 
a trial solution for the dual problem that matches the 
value Eq. To this end, we chose J^=|, Ape = Ap = 0, 

and (Ae)ij=0 except for 

(Ae)i,i = ^2, (Ae)l,8 = (Ae)8,l = ^362- 

Furthermore, 

ifJ'A)iS = (MA)i+40,i+40 = ^{l^A)i,i+iO = 

— {^J-A)i+40A = (^5)4+8, 4+8 = (/"s)i+32,i+32 = 

— (M-B)i+8,i+32 = ^(("5)4+32, i+8 = (A'c)4+16,4+16 = 

(/^c)i+24,i+24 = ^ (/"c)i+16,4+24 ~ ^ (/^c)4+24,4+16 — ^2 

for i = 9, . . . , 16. Finally, one chooses the matrices X^^ , 
A^^^ and A^^*^ . As A^^^ and A^''^ can be obtained from 
A^^'''^ by cyclic permutations, we only need to specify 
Xj^"^ . For i, j = 1, . . . , 8 we have 

iX/'^)i,j = (A^''-*)56+i,56+i = Xij, 
i^A ^)i,56+j ~ i^A ^)56+4j" = ^ij, 
(A^^"* )4+8 j+8 = (A^'''* )48+i,48+j = ^ij, 



where the nonzero elements of X and Y are given by 

XiA = 1, -^^4,4 = A^6,6 = X(i^4 = ^4^6 = 25/16, 

-^^2,3 = ^2,5 — X:^2 = -^^5,2 = ^5/4, 

Fl,! - ^4,4 - >6,6 - -1/3, 

>2,2 - ->7,7--l, 

^2,3 = Y2,5 = ^3,2 = ^3,3 = ^3,5 = —^4,6 = ^2/3, 

^5,2 = 5^5,3 = —^5,5 — —Y&A — —^8,8 — "2/3, 

>6,7 = ^7,6 = 7/80, 

^7 4 = Kt 7 = (-42 -I- \/l59559)/1200. 

A direct calculation, ideally employing a software capable 
of symbolic manipulations, now shows that these values 
determine a feasible point of the dual problem. The dual 
function for the above choice yields the value 6&2, i-e. 
the same as for the primal problem which establishes the 
optimality of the solution for the primal problem. 

APPENDIX B: FROM GHZ TO W EMPLOYING 
NON-TRACE PRESERVING PPT MAPS 

In this appendix we determine the optimal success 
probability for the transformation of a GHZ state to a 
W state under non-trace preserving CP-PPT maps. This 
problem is equivalent to the maximization of 

tr{*(pGHz)} = tT{npGHZ(^t} (Bl) 

under the constraints 

tT{npGHZ ® (1 - pw)} = 0, 

> 0, try.{r2(*)} < 1, (B2) 

This problem possesses the same symmetries (a) - (g) 
presented in section IIVI Following the same arguments 
as in section llVl most matrix elements of ft vanish. In the 
following we will present those non- vanishing matrix ele- 
ments that are sufficient to reconstruct all the remaining 
non-zero elements of the trial solution from the symme- 
tries of the problem. With bi = 0.8/3 = 262 = 463 we 
find 

f^ooiooo, 001000 = 2flooiiii, 001111 = bl, 

f^ooiooi, 001001 ~ f^ooioio, 001010 = f^ooiioo, 001100 = &2, 

f^OOlOOl, 001010 = f^OOlOOl, 001100 = —f^OOlOlO, 001100 — —b2, 

f^OOlOll, 001011 = f^OOllOl, 001101 ~ f^OOlllO, 001110 = 63, 

f^OOlOll, 001101 = —i^OOlOll, 001110 = —i^OOllOl, 001110 ~ ^3, 

f^OOOOOO, 000000 = SiloOOlll, 000111 = bl, 

f^OOOOOl, 000001 = f^OOOOlO, 000010 — f^OOOlOO, 000100 = &2, 

^^000001, 000010 = f^OOOOOl, 000100 — f^OOOOlO, 000100 = ^2, 

^^000011, 000011 — f^OOOlOl, 000101 — f^OOOllO, 000110 — bs, 

^^000011, 000101 — f^OOOOll, 000110 = f^OOOlOl, 000110 — bs, 

f^OOOOOO, 111000 — 8iloooill,llllll = — ^1, 

f^oooooi, 111001 — ^^000010, 111010 — f^oooiooaiiioo — ^2, 
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f^OOOOOl, 111010 = f^OOOOOl, 111100 = f^OOOOlO, 111100 = ^2, 
f^OOOOll, 111011 — f^OOOlOl, 111101 = f^OOOllO, 111110 = ^^3: 
f^OOOOliaillOl ~ f^OOOOll, 111110 — f^OOOlOl.lllllO — ~&3- 

Now an elementary but lengthy calculation shows that 
the chosen parameters define a feasible point of the prob- 
lem and yield a success probability of tj:{flpGHZ (E) i-} — 
0.8. 

To prove the optimality of this result we now consider 
the dual problem. The Lagrange function for the mini- 
mization problem in Eq. (|B1|I is given by 

i(fi,AA, Ab,Ac, Ap,z^) = tr{ i^pGHZ «) (1 - /Ow))} 
+tr{Ap} - tv{n{{pGHZ - Ap) ® 1 -h A^^ + A^« -I- A^^^}, 

where Ap, A^'^, A^^, A^*^ > 0. The Lagrange function has 
to be minimized over all 57 > 0. This is feasible only if 

PGHZ«)l+A^'*+A^^+A^'=-Ap(g)l-i/pG//z®(l-pvv) < 

(B3) 

in which case we obtain the dual function 

g(AA, Ab, Ac, Ap,i/) = -tr{Ap}. (B4) 

Maximizing this function under the constraints 
A^,Ab,Ac,Ap > and Eq. (jB3p yields upper bounds 
on the success probabilities of the primal problem. The 
following trial solution yields —tr{Xp} — —0.8 satisfying 
all the constraints and matching the value of the primal 
optimum thereby proving its optimality. For simplicity 
we only give the non-zero matrix elements 

Apia = \8,8 = ~\i,8 = ~Ap8,i = 0.4, 

Aai,1 = —^AlA = — A^i_6 = ^AiA = Ayi6,6 — ^A8,S 
— Aa 4,6 = 4Aa5,5 = —2Aa53 = 0.8/3, 

Aa 57,57 = —Aa 57,60 = ~Aa 57,62 — Aa 60,60 = A^ 62,62 
= Aa64,64 = Aa60,62 — 4Aa61,61 = ^ 2Aa61,64 

= 0.8/3, 
ly = 1.8. 

The elements of As and Ac are obtained from Xa by 
cyclic permutation of the parties A, B and C so that 
for example Xaz,^ = ^B2,2- Direct calculation no shows 
that this trial solution is feasible for the dual problem 
and yields the value g — —0.8 which is identical to that 
obtained from the trial solution for the primal problem. 
This completes the proof of optimality. 

APPENDIX C: FROM W TO GHZ EMPLOYING 
NON-TRACE PRESERVING PPT MAPS 

The optimization of the success probability for the 
transformation from W to GHZ proceed along very sim- 
ilar lines as those given in the previous appendix. Math- 
ematically the problem is formulated as 

tr{*(piv)} = tvinpw^t} (CI) 



under the constraints 

tr{npw ® (1 - PGHz)} = 0, 
try'{0} < 1, n^'' > 0, 

Symmetries analogous to those presented in the previ- 
ous sections hold. Following the arguments analogous to 
those in section Hvl most matrix elements of fl vanish. In 
the following we will present those non- vanishing matrix 
elements that are sufficient to reconstruct all the remain- 
ing non-zero elements of the trial solution from the sym- 
metries of the problem. With bi = 3^2/4 = 363 = 1/6 

^^000001, 000001 = f^iiiooi, 111001 = 11/90, 

^^001001, 001001 = 4iloioooi, 010001 = 4riiooooi, 100001 = ^2, 

^^001001, 010001 = f^ooiooi, 100001 = —251oioooi, 100001 = ^^2/2, 

^^011001, 011001 = f^ioiooi, 101001 ~ 2iliioooi, 110001 ~ bi, 

^^011001, 101001 = —2rioiiooi, 110001 ~ —2riioiooi, 110001 = ^3, 

^^000000, 000000 = f^iiiooo, 111000 = ^2/2, 

^^001000, 001000 = f^OlOOOO, 010000 = f^lOOOOO, 100000 = ^3, 
^^001000, 010000 = f^OOlOOO, 100000 = ^^010000, 100000 = ^3, 

^^011000, 011000 = f^ioiooo, 101000 = ^^110000, 110000 = ^3/2, 
^^011000, 101000 = f^oiiooo, 110000 = f^ioiooo, 110000 = ^3/2, 
f^oooooo, 000111 = —7riiiiooo, 111111/10 = —7/90, 

f^OOlOOO, 001111 = f^OlOOOO, 010111 = f^lOOOOO, 100111 = ^3, 
f^OOlOOO, 010111 = f^OOlOOO, 100111 = f^OlOOOO, 100111 = ^3, 

f^oiiooo, 011111 = f^ioiooo, 101111 = f^iioooo, 110111 = ^3/2, 
f^oiiooo, 101111 = f^oiiooo, 110111 = f^ioiooo, 110111 = ^3/2. 

With this trial solution we find triripw (E Pghz} = ^■ 

To prove the optimality of this result we now consider 
the dual problem. The Lagrange function for the mini- 
mization problem in Eq. IjCip is given by 

L(ri, Aa, As, Ac, Ap, I/) = -trAp - tr ^ Af'®^' 

i=A,B,C 

-tr{n{{pw - Ap) 1 - i^pw (1 - Pghz)},{C2) 

where Ap, Ayi, As, Ac > 0. This Lagrange function has to 
be minimized over all H > which is feasible only if 

{pw - Ap) (g) 1 - i^pw (1 - pghz) + Yj - 

i=A,B,C 

(C3) 

in which case we obtain the dual function 

g(AA, As, Ac,Ap,j/) = -tr{Ap}. (C4) 

Now we need to maximize this function under the con- 
straints Ayi,AB,Ac,Ap > and Eq. HC3|I . Each trial 
solution gives an upper bound on the success probability 
of the primal problem. It turns out that we can approach 
the — tr{Ap} = — | arbitrarily closely. 
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We begin by determining all non-zero matrix elements 
of Aj, in terms of Ap2,2 so that 

^p3,3 = ^p5,5 — Ap2,2, 

Ap2,3 = Ap2,5 = Ap3_5 ~ — Ap2,2/2. 

Furthermore, we completely determine the matrices 
A^, As and Ac. To this end we give all the nonzero val- 
ues of A^ as the other matrices are uniquely determined 
through cyclic permutations from A^. 



Ayl 18,23 
Ayl 18,39 



1 

5 " 
3 

10 



A 



A 34,39, 



34,23 



and 



Aai7,17 = A^33, 33 = —Aa 17,33 = —Aa33, 17 = 0.1/9, 

Aa18,18 = Aa 34,34 = 0.3; Aa 18,34 = Aa 34,18 = 0.2, 

4A^19^ig — A^35^35 — 2A^ ig 35 = 2A^35_i9 = 0.4/9, 

Aa20,20 = 4Aa 36,36 = 2Aa 20,36 = 2A^ 36,20 = 0.4/9, 

Aa2i,2i ~ 4A^ 37,37 = 2Aa 21,37 = 2A^ 37 21 = 0.4/9, 

4Ayi 32,32 — A/1 38,38 — 2Ayi 32,38 — 2A^ 38,32 ~ 0.4/9, 
Aa 23,23 = '^A 39,39 = 0.3; A^ 23,39 = 39,23 = 0.2, 
Aa24,24 = Aa 40,40 = —Aa 24,40 — —Aa 40,24 = 0.1/9. 

The elements of As and Ac are obtained from Xa by 
cyclic permutation of the parties A, B and C so that for 
example Xa5,5 — As 2.2- A direct calculation shows that 
the constraints Xa, As, Ac, Ap > are satisfied with these 
choices. Now we need to verify whether the constraint 

{pw -Xp)(g)l- vpw «) (1 - PGHZ) + - 

i=A,B,C 

(C5) 

can be verified as well. Note that we still have the free 
parameters Ap2,2 and v. A lengthy computation (prefer- 
ably employing Mathematica) shows that the left hand 
side of the constraint has 6 distinct nonzero eigenvalues, 
namely 



M2 = ^(13 - 135Ap2,2), 



Hi = 2-v, 

P3 = -^(-2 - 45Ap2,2), Pi = ^(4 - 45Ap2,2), 
P± = ^[47-30i.-45Ap2,2± 



1569 - 2220// + 3330Ap2,2 + (45Ap2,2 - 30j/)2 ]. 



Clearly, for > 2 and Ap2,2 > the first 4 eigenvalues 
are non-positive. Now we can verify by direct inspection 
that for any choice of Ap2,2 > 1/9 there is a choice of 



1/ > 2 such that the two eigenvalues /z± are negative so 
that also the constraint Eq. (|C5() is satisfied. There- 
fore, for any value of — tr{Ap} < — | we can satisfy the 
constraints. This shows that the primal problem which 
achieves a success probability p — 1/3 is optimal. 



APPENDIX D: SINGLE COPY DISTILLATION 
FROM HIGH RANK MIXED STATES 

In this appendix, we prove that PPT-operations can- 
not distill any pure entangled states from a single copy of 
p on C''(8'C'^ when rank{p) > (P—2. To this end, it suffices 
to show that the success probability p{p P^) under 
PPT-operations (5') in the trace non-preserving scheme 
is strictly zero, where p&T-l{V) and &H{V'). Since 
both p^P^ and (1— P^) is invariant under the local 
unitary transformation oi 1 ^ 1 ^ U ^ U* , it suffices to 
consider Q invariant under these local operations, i.e. 



A(g)P;^ + B< 



1 



(Dl) 



with A and B being matrices on 7i(V^). The success 
probability is then 



Pip ^ P+) = tT{np ® P+} - tr{^p}, 
and constraints for Q are 



(D2) 



tr{fi/9 ( 
A>0, 
1 



(t 



B>0, 



-B'^>A'^> 



= tr{Bp} = 0, 
1>A + B, 



1 



Since B>0 and tr_Bp = 0, the support space of B must be 
contained in the kernel space of p, and hence rank(i?) < 
2 when rank(p) > tP — 2. On the other hand, B^'^ > 
must hold from -77^ B^^ > --rj^B^-^, and B must 
be a separable state (leaving out normalization) since 
rank(_B) <d Therefore, by using appropriate local 

basis, B can be written as 



i? = y|ll)(ll| + z|e/)(e/|, 
where y and z are non-negative values and 



(D3) 



\ef) = (cosu|l) + sinu|2)) (g) (cosi;|l) -|- sinu|2)) (D4) 

is a product vector. In this choice of local basis, B^'^ =B. 
Let P be the projector on the support space of B^^^ and 
Q = I-P. The condition of jt^B^^ >A^^>--^B^^ 
implies that ±QA^^Q>0, and hence QA^^Q = must 
hold. Furthermore, A^"^ + -jq-^B^^ must be a positive 
operator, for which Q{A^^ + ■^p^B^^)Q = also holds. 
Therefore, support space of A^ + -jr^i B^ ^ must be P, 
and hence the support space of A^^^ must be contained 
in the support space of B^^. As a result, rank(A'"'*) < 
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rank(_B'"-*) < 2. Furthermore, A^^ must be written in 
the form of 

A^- = r|ll)(ll| + ,s|ll)(e/| + .s*|e/)(ll| + <|e/)(e/|, 

and A is then given by 

A = r|ll>(ll| + ,s|el>(l/| + s*|l/>(el| + t\ef){ef\. 

Therefore, A must be essentially two-qubit state (leaving 
out normalization) since A> must hold. If the two- 
qubit state A is entangled, rank(A^^) must be 4 |4^.l49|. 
which contradicts that rank(A^^)<2. Therefore, A and 
A^^ must be written in a separable form. 



In the case where sin u sin w 0, the support space of 
A^^, which is spanned by |11) and |e/), contains only 
two product vectors (|11) and \ef) itself) and hence 
A^^ must be written as 

^r|ll)(ll|+t|e/)(e/|, (D5) 

and A — A^^. As a result, the support space of A is 
contained in the support space of B, and hence p{p — > 
P^) = trAp= astri?p = 0. In the case where sin m sin 11 = 
0, |e) = |1) or I/) = |1) holds. As a result, A is spanned 
by {|11), |1/)} (or {|11), |el)}) which is a kernel of p, and 
hence p{p^ P^)=trAp~0. 
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